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RESEARCH MEMORANDUM 


EXPERIMENTAL AND ANALYTICAL INVESTIGATION OF THE TRANSONIC 

AND SUPERSONIC DIVERGENCE CHARACTERISTICS OF A 

DELTA-PLAN-FORM ALL-MOVABLE CONTROL* 

By A. Gerald Rainey, Perry W. Hanson, and 
Dennis J. Martin 


SUMMARY 


The static aeroelastic divergence characteristics of a delta-plan- 
form model of the canard control surface of a proposed air-to-ground 
missile have been studied both analytically and experimentally in the 
Mach number range from 0.6 to 3 . 0 . The experiments indicated that 
divergence occurred at a nearly constant value of dynamic pressure at 
Mach numbers up to 1.2. At higher Mach numbers somewhat higher values 
of dynamic pressure were required to produce divergence. The analysis 
and the experiment indicate that the camber stiffness of the control 
surface and the stiffness of the control actuator are both important in 
divergence of surfaces of this type. 


INTRODUCTION 


The increased usage of low-aspect-ratio canard surfaces for stabil- 
ity and control of missiles has led to considerable interest in the 
aeroelastic characteristics of such surfaces. In several instances, 
missile failures have occurred which were believed to be due to static 
aeroelastic divergence of surfaces of this type. In most cases, a 
relatively simple solution to the problem has been found such as 
stiffening the surface in the chordwise direction or by altering the 
geometry of the control. Investigations of this type are usually 
of an ad hoc nature and the results may not be generally available. 
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Divergence encountered by surfaces of this type differs somewhat 
from the classical torsional divergence in that camber deformations 
seem to play a dominant role. This new class of divergence problems 
which is associated with thin low-aspect-ratio surfaces has received 
some analytical study (refs. 1 and 2 ). 

In a recent design of an air-to-ground missile, preliminary studies 
showed that the canard control surface might be subject to divergence 
within the flight boundary of the missile. Consequently, a series of 
models were constructed and have been tested in the Langley 2-foot tran- 
sonic flutter tunnel in the Mach number range from 0.6 to 1.2 and in 
the Langley 9 - by l 8 -inch supersonic flutter tunnel over the Mach num- 
ber range from 1.64 to 3.0. In addition to a simple proof test of the 
missile configuration as designed, the opportunity was taken to investi- 
gate the effects of variations of stiffness and location of the pitch 
axis. An analytical treatment of the divergence of this type of control 
has been developed. The structure has been treated as a beam with its 
span alined with the airstream. Two types of aerodynamic forces are 
considered, one based on very- low-aspect-ratio theory and the other 
based on piston theory. The experimental results are compared with the 
results of this analysis to aid in making the investigation of more 
general interest. 


SYMBOLS 


A - * slope influence coefficient for panel, pitch spring being 

considered infinitely stiff (slope at position i due to 
unit load at position j), radians/lb 

a. . slope influence coefficient for pitch spring, panel being 

considered infinitely stiff, radians/lb 

a speed of sound, ft/sec 

b model semichord measured parallel to the root chord at three- 

quarter span, ft 

B. . slope influence coefficient for panel-spring combination 

1J ( B ij = a ij + A ij)' radians/lb 

c distance from intersection of leading edge and root chord 

to trailing edge, ft 
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measured deflection of control surface at a point 27. 1 percent 
of the root chord rearward of the leading edge and at 50 per- 
cent of the local span outboard of the root chord due to a 
unit load at that point with pitch spring stiffness assumed 
infinite, ft /lb 

modulus of elasticity of panel, lb/sq ft 
modulus of elasticity of air, lb/sq ft 


effective value of modulus of elasticity of control, lb/sq ft 

p 

panel bending stiffness with respect to pitch axis, lb-ft^ 

measured deflection of infinitely stiff control surface at a 
point 27-1 percent of the root chord rearward of the leading 
edge and at 50 percent of the local span outboard of the 
root chord due to a unit load at that point acting against 
the pitch spring stiffness only, ft/lb 


effective stiffness of panel-spring combination, 
lb/ft 


1 

d + h* 


length of trailing edge, ft 
bending moment, ft-lb 


mass of panel, slugs 

K e 

effective mass of panel, — , slugs 


mass of air contained in the cone whose base diameter is equal 
to the root chord and whose height is equal to the span, slugs 

M Mach number 

p static pressure, lb/sq ft 

P aerodynamic load, lb 

q dynamic pressure, lb/sq ft 

s one-half the distance from leading edge to root chord meas- 

ured parallel to pitch axis at chordwise station x, ft 
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thickness of panel, ft 
stream velocity, ft/sec 

component of stream velocity normal to control surface,- ft/sec 

chordwise distance measured from and perpendicular to pitch 
axis (for analytical purposes, pitch axis is assumed perpen- 
dicular to line bisecting angle formed by leading edge and^ 
root chord), ft 

chordwise station where deflection is measured due to load 

at x*, ft 
J 

distance of panel elastic axis from leading edge at root 
chord, ft 

chordwise station where load is placed, ft 

vertical displacement, ft 
ratio of specific heats 
Mp 

mass ratio, — 

m a 

density of air, slugs/cu ft 

natural frequency of vibration, radians/sec 

spring constant of pitch spring, ft-lb/radian 




distance from body center line to control leading edge, ft 
= 1 if i = j 
= 0 if i i j 

M differentiating matrix 

Subscripts: 

L refers to lower surface 
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U refers to upper surface 

00 refers to conditions far removed from control surface 


APPARATUS AND TESTS 


Description of Wind Tunnels 

The tests were conducted in the Langley 2-foot transonic flutter 
tunnel for the Mach number range from 0.6 to 1.2 and in the Langley 
9- by 18- inch supersonic flutter tunnel for the Mach number range 
from 1.64 to 3.0. 

The Langley 2-foot transonic flutter tunnel is a slotted- throat 
single-return wind tunnel equipped to use either air or Freon- 12 as a 
test medium. All of the present tests were made with Freon-12. The 
tunnel is of the continuous-operation type, powered by a motor-driven 
fan. Both test -section Mach number and density are continuously 
controllable . 

The Langley 9- by l8-inch supersonic flutter tunnel is a fixed- 
nozzle blowdown- type wind tunnel exhausting into a vacuum sphere. The 
nozzle configurations used in this investigation gave Mach n umb ers of 1.64, 
2.0, 2.55, and 3*0. At each Mach number the test-section density varies 
continuously to a controlled maximum. 


Description of Models 

The l/9-scale models simulated the delta plan form of the canard 
all-movable control surfaces of an air-to-ground missile. They were cut 
from 2024-T aluminum sheet stock, the thickness of a given model being 
constant over the plan form except for the beveled leading and trailing 
edges. The geometry of the models and model-mount fairings is shown in 
figure 1. The portion of the mount fairings forward of the trailing edge 
simulated the contour of the missile. 

The masses and thicknesses of the control- surf ace models, identi- 
fied by numbers 31 to 47, are presented in table I. The method of 
mounting the models for use in both the 9- by l8-inch supersonic flut- 
ter tunnel and the 2-foot transonic flutter tunnel is shown in figures 2 
and 3- The torque rod was connected to the mount frame through a tor- 
sional spring. Several torsion springs were used to cover a range of 
stiffnesses. Basic combinations of torsion springs and control- surface 
thickness produced model modes simulating the symmetrical and antisym- 
metrical modes of the prototype control surface. In addition to the 
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basic combinations several modified combinations were used to increase 
the scope of the investigation. It should be noted that although the 
physical appearance of the model mounts was different; the model root 
conditions were the same in both mounts . 

A model is shown mounted in each of the tunnels in figures 4 and 5. 
Also shown in figures 4 and 5 are the different mount fairings used in 
the two tunnels. The differences in model mount fairings are also indi- 
cated in figure 1. 

The torsional stiffnesses of the springs are presented in table II, 
along with the control- surface-panel stiffnesses and combination panel- 
spring stiffnesses. Thus, model 35-55 is control surface model 35 
mounted on spring 55. The column headed d is the measured deflection 
of the control surface at a point due to a unit load at that point with 
the pitch spring stiffness assumed infinite and the column headed h 
is the deflection of an infinitely stiff control surface at a point 
due to a unit load at that point acting against the pitch spring stiff- 
ness only. The point of reference is at 27. 1 percent of the root chord 
rearward of the leading edge and at 50 percent of the local span out- 
board of the root chord. The effective stiffness K e is a measure of 

the total stiffness of the model and is defined as — — ■— pounds per 

d + h 

foot. Also shown in table II are calculated divergence dynamic pres- 
sures obtained from an analysis to be discussed subsequently. 


Test Procedure 

Tanpiey 9- by l8-inch supersonic flutter tunnel .- The models tested 
in the Langley 9- by l8-inch supersonic flutter tunnel were all of the 
basic configuration; that is, the spring and control surface combinations 
were such that the elastic properties of the actual canard all-movable 
control were simulated, as was the location of the pitch axis ( 0.62 root 
chord). Electrical resistance wire strain gages were mounted at the 
root near the hinge line and the signal was taken to a recording oscil- 
lograph which also recorded tunnel conditions. In addition, high-speed 
motion-picture cameras recorded the behavior of the model. The proce- 
dure for making all the runs was as follows: the models were set at 

zero angle of attack and then the tunnel was evacuated to approximately 
1 in. Hg absolute. A control valve upstream of the test section was 
then opened and the density of the flow was allowed to increase at 
constant Mach number until divergence occurred. 

Langley 2-foot transonic flutter tunnel .- In addition to the basic 
configuration, several modified configurations were tested in the 
Langley 2-foot transonic flutter tunnel. Effects of variation of the 
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pitch-axis location, variations of spring and control- surf ace stiffnesses 
were investigated. In order to obtain data at various Mach numbers, the 
following procedure was used. With the tunnel set at a low density the 
velocity was increased until the desired Mach number was reached. With 
the velocity held approximately constant, the test-section density was 
slowly increased until divergence occurred. The dynamic pressure was 
then decreased rapidly by actuating a spoiler in the diffuser section 
of the tunnel. The Mach number was then decreased to a point well below 
the divergence condition. At this point the stagnation pressure was 
increased by a small amount, then the velocity was slowly increased 
until divergence occurred. This procedure was repeated for several small 
increments in stagnation pressure. For the type of boundary found for 
these models this procedure resulted in divergence points for several 
Mach numbers from the maximum obtainable in the tunnel down to some arbi- 
trary lower Mach number. 


Data Reduction 

It was necessary to test models of different stiffnesses in order 
to obtain divergence data over the desired range of Mach number within 
the range of dynamic pressure obtainable in the test facilities. This 
variation in stiffness leads to the necessity of reducing the data 
obtained for the various models to some form of dimensionless parameter 
which will provide a basis for comparison of the test results at various 
Mach numbers. Such a parameter has been developed and discussed in 
appendix A. The parameter chosen is closely related to the stiffness- 
altitude parameter which has proven useful in interpreting flutter 
results. The divergence parameter differs from the flutter parameter 
in that the frequency and mass have been replaced by a stiffness term 
in an attempt to recognize the static characteristics of divergence. 

This parameter is 

b K 
a V H A 

where b is a reference semichord taken at the 75-percent- semispan 
station, a is the speed of sound, and K e . is the deflection stiff- 
ness or the load required for a unit deflection measured at an arbi- 
trary point on the surface. For all the models tested, b is 
0.0926 foot, and is the mass of air contained in a volume of a cone 

whose base diameter is equal to the root chord and whose height is equal 
to the exposed span of the control surface. This volume is 0.0^05 cubic 
foot. 
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The results of the experiments will he discussed subsequently, along 
with results obtained from the following analysis. 


ANALYSIS 


This section is concerned with the development of divergence equa- 
tions applicable to the spring-mounted elastic control tested in the 
experimental investigation. An influence -coefficient method of analysis 
is used in which two different methods are used for representing the 
aerodynamic forces, namely, low-aspect-ratio theory (ref. 3), and piston 
theory (ref. 4). 


Structural Representation 

In order to structurally represent the surface in a manner that is 
readily amenable to analysis, the sections of the surface were considered 
to be sheared parallel to the pitch axis and the trailing edge was rota- 
ted about its midspan point so that an equivalent symmetrical plan form 
was obtained. The equivalent plan form is indicated in the following 
sketch: 



In both the low-aspect-ratio-theory and the piston- theory approaches, 
the aerodynamic loading is defined in terms of the local streamwise 
slopes and curvatures. The expressions for aerodynamic loading can be 
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combined with the slope influence coefficients of the system to obtain 
an expression for the dynamic pressure at divergence. A first step in 
the development of the divergence equations is the determination of the 
combined slope influence coefficients of the spring-mounted elastic sur- 
face. The simple beam equation will be applied in the stream direction 
to determine a slope influence coefficient array. . It is assumed that 
the influence of spanwise deformations on the structural and aerodynamic 
forces is small. The elastic influence coefficients of the surface alone 
may be obtained by assuming that C^, the pitch spring, is infinitely 

stiff; that is, slope and deflection at the pitch axis are zero. Use 
may then be made of the fundamental beam relation 


dx 


(i) 


For a concentrated load P v applied at a point at a distance x. from 

x 0 J 

the pitch axis, equation (l) becomes 


ei £t - • h 


x. - x)P 


( 2 ) 


for 


x < 


X . 
0 


Since the surfaces considered are of constant thickness, the section 
moment of inertia I may be written as 


1 - S? (*p - *) 


(3) 


Equation (2) may be integrated with the section moment of inertia repre- 
sented by equation (3) to obtain the slope at a point x^ due to a 

load at station x. 



12ncP„ 

A 
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where 

n = 1 
n = 0 
n = -1 

Appropriate boundary conditions are 


( 0 < h<x 3 ) 

h - °) 

h <x i <o ) J 



(in 


i - Vi > 0) 



( x i x j 


< 0 ) 




(5a) 


(5b) 


An elastic slope influence coefficient A^ may then be defined as, 



subject to the conditions of equations (5)* 

For the present analysis the control surface is divided into ten 
sections of equal increments along x and the control points are located 
at the middle of each section. The ten-point slope influence coefficient 

matrix [a] calculated from equation (6) and representing the control 
surfaces is presented in appendix B. The notation LI represents a 
square matrix. 


The slope influence coefficient a^ due to a spring in the pitch 
degree of freedom is 


x i 

= c,, 


(7) 


The matrix [a] representing the pitch springs is also presented 
in appendix B. 
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The combined slope influence coefficients due to the elastic con- 
trol surface and the spring restraint in the pitch degree of freedom 
are additive, 


B ij ” A ij + a ij 


(8) 


The matrix equation 



(9) 


gives the slope <j^j> 
loads ^ P | , where j- 


in the streamwise direction for any system of 
denotes a column matrix. If the aerodynamic 


loads can he expressed in terms of the dynamic pressure and slope, then 
substitution of the aerodynamic loads into equation ( 9 ) results in the 
divergence equations which may be iterated to obtain the critical dynamic 
pressure. Two methods of representing the aerodynamic loads will be 
used, namely, low-aspect-ratio theory and piston theory. The following 
section presents the development of the aerodynamic loads into a form 
that can be used in equation ( 9 ) to obtain the divergence equation. 


Divergence Equations 

Low-aspect-ratio theory .- The aerodynamic loads are first obtained 
from very-low-aspect-ratio theory (ref. 3)* This theory assumes that 
the flow field within a planar strip perpendicular to the flow direc- 
tion is two dimensional and that the changes in the flow direction are 
small. The complete expression for the aerodynamic load on a section 
of dimension 2s normal to the flow and Ax parallel to the flow may 
be written as 

P = -rtp(Ax)s 2 ^ + zV ~ + V 2 j - 2n(Ax)pVs tan e(z + V 

( 10 ) 

where 0 is the angle that the leading edge is inclined to the free 
stream. The effects of the central body on the aerodynamic forces as 
given by low-aspect-ratio theory are not known; however, they are assumed 
to be small. The time derivatives for the divergence case vanish, and 
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equation (10) when applied over the entire control surface may be 
written in matrix notation as 


(P} = -2*(A*h | js 2 g| + 2 tan 8 js g| 


If 


( 11 ) 


5 ij - 1 

(i = 

s ij = 0 

(i t 


(12) 


the matrix for the aerodynamic loads becomes 


jp| =2n(Ax)q<j 


6. .s / 
ij J 


d^z 

dx 2 


> + 2 


tan 0 [b. .s .1 I — | 
L Jj I dx | 


A differentiating matrix |Vj may be determined such that 


siJ 

dx S 


= M 1^1 

.2 | * — ' ^ dx J 


(13) 


(14) 


A sample matrix jV] for the ten-point analysis used in this paper is 
given in appendix B. If the differentiating matrix []d] of equation (l4) 
is used in equation (13)> the expression for the aerodynamic loads 
becomes, 


P = -2rt(Ax)q 


5 « S J 2 ] M + 2 tan e 


/ dz \ 


(15) 


dz 


The square matrix premultiplying < — > is a function of geometry only, 

L dx J - “ 

and, if it . is denoted by [C J , the aerodynamic loads are given by 

{p} = -2«(toh[c] |gj> (16) 
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If the aerodynamic loads given by equation (l6) are substituted into 
the combined slope influence equation (eq. (9)), the equation governing 
the slopes under aerodynamic loadings are 

{g} =- 2 ^hrB][c]{g} (17) 


Equation (17) expresses the conditions for which the aerodynamic forces 
are equal to the structural restoring forces. Equation (17) is thus the 
divergence equation and may be iterated to obtain the dominant root which 
yields the dynamic pressure at divergence. The values of q thus obtained 
for each case are given in table II. It may be mentioned that the pro- 
duct [b] [c] for stiff control surfaces and weak pitch springs produced 

an ill-conditioned matrix which was divergent under normal iteration pro- 
cedures. Averaging successive iterations proved to be adequate to force 
convergence to the dominant mode in the cases treated. 

Piston theory .- A second method of representing the aerodynamic 
forces for the supersonic case was also used and involved the use of 
piston theory (ref. 4). Piston theory is an application of the "local" 
wave equation and may be obtained from potential-flow theory if the 
Mach number is allowed to take on large values . The pressure coeffi- 
cient may be written as 


p - p = p a 
* -^00 -^00 00 


w + ( 7 + l\ /w \2 / y + l \ / w 

a °o \ 4 A 8-00 ' \ 12 A a oo/ 


+ • . 


(18) 


The load on a section of the upper surface, which is 2s wide and Ax 
long, becomes 


Py = -2(Ax)(2s)q 


i© + (HA# + AHA# 


(19) 


The surface is of constant thickness and the load on a section of the 
lower surface is 

P L - -2(Ax)(2s)q jj. $) + . M 5 (^)Gf + . . .] 


( 20 ) 


Recognizing that dz/dx is equal to w/v, the total load P - 
becomes 
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P = -4(Ax)(2s)q 


— — + W* 
M dx 




( 21 ) 


Only the first term of equation (2l) will be used in the present analy- 
sis. The system of equations representing the loads on the control sur- 
face is 



( 22 ) 


Equation (13) is the corresponding equation derived from low-aspect- 
ratio theory. The square matrix premultiplying in elation (22) 

is also a function of geometry only, and, if it is denoted by [c] , the 
aerodynamic loads are 


{Pj = - i (Ax)q[c] 


(23) 


Substituting the aerodynamic loads given by equation (23) into equa- 
tion (9) gives the divergence equation for the analysis based on piston 
theory 


{ 


dz 

dx 


= - 4 (ta)q[B] 0 {g} 


(24) 


Equation (24) may be iterated to obtain the critical values of q. 


RESULTS AND DISCUSSION 


The basic model configuration with springs simulating both the . 
symmetric and the antisymmetric stiffness of the prototype missile has 
been tested in the two wind tunnels in the Mach number range, from 
about M = 0.6 to 3.0. Additional tests have been made in the. tran- 
sonic tunnel to study the effects of stiffness of the control rotation 
springs and of the control surface. Additional studies were made of 
the effects of location of the pitch axis. 
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General Characteristics of the Divergence Encountered 

Classically, divergence has teen treated as an aeroelas tic phenom- 
enon associated with torsional deformations. This type of divergence 
has' been defined as a static instability of an airfoil in torsion which 
occurs when the torsional rigidity of the structure is exceeded by 
aerodynamic twisting moments (ref. 5 ) • The type of divergence encoun- 
tered in the present investigation seems to fit this same definition 
except that the role of torsional deformations has been replaced by 
camber deformations superimposed on a rotation of the control about 
its pitch axis. The type of motion involved is shown in figure 6 
which is composed of enlargements from a high-speed motion picture. 

As the deflections become large it can be seen that the surface has 
large curvature ahead of the pitch axis and a decided slope at the 
pitch axis. As a matter of interest, deflections were me asured on 
several of the enlargements and are compared to the calculated deflec- 
tion shape in figure 7. The agreement between the measured and cal- 
culated deflection shapes is good. 

The type of motion involved in divergence of these models is quite 
violent in the sense that very large deflections are reached in a very 
short period of time as indicated by the enlargements of the high-speed 
motion picture shown in figure 6. At subsonic and transonic speeds 
only a few of the models acquired a permanent set during divergence, 
presumably because of a stalling effect at high angles of attack. At 
supersonic speeds , all the models were permanently damaged in divergence . 
A representative selection of these damaged models is shown in figure 8. 
Although the models did not always suffer damage at the lower Mach num- 
bers the control deflections during divergence were probably sufficiently 
large to cause very violent maneuvers of the missile with subsequent 
structural damage . 


Basic Configuration 

The data obtained for the basic configuration have been reduced 
to a nondimens ional stiffness-altitude parameter which is discussed in 
appendix A. The values of this parameter represent a stability boundary 
for static aeroelastic divergence and are shown as a function of Mach 
number in figures 9 and 10. In a figure of this type, constant altitude 
operation of a given configuration would be represented by a horizontal 
line at a value of the parameter determined by the stiffness of the 
control and the altitude. Radial lines through the origin represent 
lines of constant dynamic pressure. 
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The tendency toward a decrease in slope of the boundary with 
increasing Mach n umb er indicates that somewhat higher dynamic pressures 
would be required to produce divergence at higher Mach numbers than at 
lower Mach numbers. 

The analysis of static aeroelastic divergence using very-low-aspect- 
ratio aerodynamic theory yields a single value of dynamic pressure 
required to produce divergence regardless of the Mach number . If pis- 
ton theory is used, the analysis indicates that the dynamic pressure 
at divergence increases directly with Mach number. The calculated 
results obtained from both types of aerodynamic theory are shown on fig- 
ures 9 and 10. In the Mach number range from about 0.6 to 1.2 (where 
piston theory is not applicable), the agreement between the experiment 
and calculations based on low-aspect-ratio theory is considered to be 
excellent. At higher supersonic Mach numbers, the experimental results- 
fall about one-half the distance between the calculated results obtained 
for the two types of aerodynamic theory. 


Effects of Variations in Stiffness 

In obtaining data over the desired range of Mach number in the two 
facilities it was necessary to use models of varying stiffness. An 
impression of the effects of stiffness can be obtained by examination 
of figures 9 and 10 and observing the degree to which a single curve 
can be fitted to the data for models of various stiffness levels. The 
fitting of a single straight line to the data implies that the dynamic 
pressure required for divergence is essentially directly proportional 
to the stiffness. This seems to be true for cases where the contribu- 
tions of the control surface and the pitch spring to the total stiffness 
remain in about the same proportion. When the relative contributions 
of the two sources of stiffness are varied, this direct relationship 
between dynamic pressure and stiffness cannot be expected to apply. 

This feature is illustrated in figures 11 and 12, where the variation 
of the dynamic pressure required for divergence with stiffness is shown 
for two methods of varying the overall stiffness of the model. The 
first method (fig. 11) was to test the same control surface mounted on 
different springs simulating a variation in control actuator stiffness. 
The second method (fig. 12) was to test control surfaces of varying 
stiffness mounted on the same spring. 

The data agree very well with the calculated values and indicate 
that the stiffness of the surface and the stiffness of the control actu- 
ator are both important in determining the divergence characteristics 
of controls of this type. 
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Effects of Variations of the Pitch Axis 

It has long been recognized that the relative location of the aero 
dynamic center of pressure and the elastic axis is important in aero- 
elastic problems. In the present investigation it was believed that 
the camber deformations of the surface were producing a more forward 
location of the center of pressure than would be the case for a more 
rigid surface and, consequently, it was considered desirable to deter- 
mine the effects of moving the elastic axis or the pitch axis forward. 
For a model which simulated the prototype stiffness of both the surface 
and the actuator, it was found that moving the pitch axis forward from 

0 . 62 c to 0 . 58 c increased the dynamic pressure at divergence by about 
55 percent. Similar tests with a much lower simulated actuator stiff- 
ness indicated about an 80 -percent increase in dynamic pressure for the 
same change in axis location. When the same control surface was tested 
with the axis at midchord and with zero actuator stiffness (free 
floating), the dynamic pressure at divergence was increased by about 
20 percent indicating the strong influence of the location of the pitch 
axis. 


CONCLUSIONS 


Divergence studies of a delta-plan-form all-movable control in the 
Mach number range from 0.6 to 5-0 indicate the following conclusions: 

1. At Mach numbers from 0.6 to 1.2 divergence occurs at an almost 
constant value of dynamic pressure. At higher supersonic speeds up to 
a Mach number of 5*0, divergence occurs at somewhat higher values of 
dynamic pressure. 

2. Analytical results based on very-low-aspect-ratio aerodynamic 
theory gave very good agreement with the experimental results in the 
Mach number range from 0.6 to 1.2. At higher Mach numbers the experi- 
mental results fell about one-half the distance between two sets of 
calculated results based on low-aspect-ratio theory and piston theory. 

5 . The analysis and the experiment indicate that the stiffness 
of the control surface and the stiffness of the control actuator are 
both important in divergence of controls of this type. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 1^, 1958. 
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APPENDIX A 

DERIVATION OF A PARAMETER FOR PRESENTATION 
OF EXPERIMENTAL DIVERGENCE DATA 


In the study of dynamic aeroelastic phenomena or flutter it has 
been found that a convenient grouping of parameters called the stiffness- 
altitude parameter has been very useful in interpreting experimental 
flutter data obtained for a variety of stiffnesses over a range of alti- 
tude and Mach number. This flutter parameter consists of the product 
of a reduced frequency based on a representative chord, natural frequency, 
and the speed of sound times the square root of a mass ratio which is 
usually taken as the ratio of the mass of the surface to the mass of a 
specified volume of air surrounding the surface. This flutter parameter 

can be written as 

If it is reasoned that static aeroelastic phenomena, in particular 
divergence, do not depend on inertia forces, then it seems logical that 
some other combination of parameters might be more useful in interpreting 
divergence data. If the divergence model can be represented by a con- 
centrated mass which yields the frequency co when attached to a spring 
whose spring constant is K e , the flutter parameter might be redefined 
as 



This new parameter would seem to be more appropriate for diver- 
gence studies since it is not based on dynamic properties of the model 
but does include the stiffness of the surface. However, the new param- 
eter is somewhat unsatisfactory because the significance of the individ- 
ual parts of the parameter is not obvious. As a matter of interest 
the parameter can be reduced further to 
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p 

where the product pb l is proportional to the mass of a particular 
volume of air surro und in g the surface. The speed of sound can be elimi- 


nated by the relationship 



E. is the modulus of elas- 
A 


ticity of the medium. Recognizing that the spring constant K g is 

proportional to an effective value of the modulus of elasticity of the 
material E^ e , the parameter becomes 



where C is a constant for a given configuration depending only on the 
geometry of the configuration. Thus, it is seen that the divergence 
parameter is, essentially, the ratio of the model stiffness to the air 
stiffness which would seem to be a very significant parameter. 


The divergence boundary defined by the dimensionless stiffness- 
altitude parameter can be converted easily to a boundary in terms of 
dynamic pressure and Mach number for a particular configuration. At 
each point on the boundary the dynamic pressure at divergence can be 
found from the following relation: 



where — is the specified volume of the medium surrounding the 
surface , 
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APPENDIX B 


SAMPLE DIVERGENCE CALCULATION 


Presented in this appendix is a sample calculation of the dynamic 
pressure at divergence. Low-aspect-ratio theory will he used. The 
control surface was represented by the following parameters. The dimen- 
sions are given for the equivalent control surface after the sections 
were sheared parallel to the pitch axis and adjusted to obtain a sym- 
metrical control surface. 

E = 10,000,000 lb/sq in. = 1,4k), 000, 000 lb/sq ft 
c = 8.55 in. = 0.7125 ft 

l = 3.60 in. = 0.50 ft 
Xp = 5.45 in. = 0.454 ft 

(Ax) = 0.855 in. = 0.071 ft 
t = 0.052 in. = 0.00267 ft 

l/Cp = 0.0545 radian/in- lb = 0.4l4 radian/ft-lb 
9 = 50° = 0.524 radian 

The influence-coefficient matrix for all the control surfaces is 
calculated from equation (6). 

3-937 2.320 1.323 .674 .271 .058 

3.554 2.320 1.323 .674 .271 .058 

2.922 2.122 1.323 .674 .271 .058 

2.216 1.702 1.188 .674 .271 .058 

r A -I _ 12c 1.475 1.174 .873 .572 .271 .058 

L iJJ E 7 t 5 -7U -580 .450 .319 .189 .058 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 


The slope-influence-coefficient matrix for the pitch degree of freedom 
is given as 


’5.019 4.170 3.322 2.473 1.625 .776 -.072 -.921 -1.769 -2.618] 

5.019 4.170 3.322 2.475 1.625 -776 -.072 -.921 -1.769 -2.618 

5.019 4.170 3.322 2.473 1.625 .776 -.072 -.921 -1.769 -2.618 

5.019 4.170 3.322 2.473 1.625 .776 -.072 -.921 -1.769 -2.618 

r ft -1 _ 1_ 5-019 4.170 3-322 2.473 1.625 .776 -.072 -.921 -1.769 -2.618 

L iJJ c p 5.019 4.170 3.322 2.473 1.625 .776 -.072 -.921 -1.769 -2.616 

5.019 4.170 3.322 2.473 1.625 -776 -.072 -.921 -1.769 -2.618 

5.019 4.170 3.322 2.473 1.625 -776 -.072 -.921 -1.769 -2.618 

5.019 4.170 3.322 2.473 1.625 .776 -.072 -.921 -1.769 -2.618 

5.019 4.170 3.322 2.473 1.625 .776 -.072 -.921 -1.769 -2.618 
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The combined slope-influence-coefficient matrix = A^j + a^ is 

given by adding equations (Bl) and (B2) 


N ■ 


.513 

•344 

.229 

.144 

.080 

.032 

-.003 

-.032 

-.061 

-.090 

.480 

.541* 

.229 

.144 

.080 

.032 

-.003 

-.032 

-.061 

-.090 

.1*25 

.327 

.229 

.144 

.080 

.032 

-•003 

-.032 

-.061 

-.090 

.361* 

.291 

.217 

.144 

.080 

.032 

-.003 

-.032 

-.061 

-.090 

.300 

.21*5 

.190 

.135 

.080 

.032 

-.003 

-.032 

-.061 

-.090 

.235 

. 19 !* 

• 153 

.113 

.072 

.032 

-.003 

-.032 

-.061 

-.090 

• 173 

.144 

.115 

.085 

.056 

.027 

-.004 

-.033 

-.063 

-.093 

• 173 

.144 

.115 

.085 

.056 

.027 

-.004 

-•038 

-.079 

-.119 

• 173 

.144 

.115 

.085 

.056 

.027 

-.oo4 

-.038 

-.084 

-.153 

• 173 

.144 

.115 

.085 

.056 

.027 

-.oo4 

-.038 

-.048 

-.138 


(B3) 


A differentiating matrix is obtained by applying the 5-point interpolation 
equations given on page 97 of reference 6. 


M- 


12(Ax) 


25 

-48 

36 

-16 

3 

0 

0 

0 

0 

0 

3 

10 

-18 

6 

-1 

0 

0 

0 

0 

0 

-1 

8 

0 

-8 

1 

0 

0 

0 

0 

0 

0 

-1 

8 

0 

-8 

1 

0 

0 

0 

0 

0 

0 

-1 

8 

0 

-8 

1 

0 

0 

0 

0 

0 

0 

-1 

8 

0 

-8 

1 

0 

0 

0 

0 

0 

0 

-1 

8 

0 

-8 

1 

0 

0 

0 

0 

0 

0 

-1 

8 

0 

-8 

1 

0 

0 

0 

0 

0 

1 

-6 

-18 

-10 

-3 

0 

0 

0 

0 

0 

-3 

16 

-36 

48 

-25 


(b4) 


The matrix [b^sjj is obtained from the geometry of the control and 
is expressed in inches as follows 


6, .s. = 
ij J 


3-6 

40 


0 

3 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

5 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

7 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

9 


0 

0 

0 

0 

0 


0 11 
0 0 


0 

0 

■0 


0 

0 

0 


0 

0 

0 

0 

0 

0 

13 

0 

0, 

0 


0 

0 

0 

0 

0 

0 

0 

15 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

17 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

19 


(B 5 ) 


also 
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3 ij s J 


■fctf 


0 

9 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

25 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

49 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

81 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

121 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

169 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

225 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

289 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

361 


(B6) 


The matrix [c] is obtained from equations (15) and (l6) 


c = 




[c] 


sj 2 ] H 

+ 2 

tan 0 




.0 84 

.038 

-.029 

.013 

-.002 

0 

0 

0 

0 

0 

-.021 

.240 

.129 

-.043 

.007 

0 

0 

0 

0 

0 

.020 

-.159 

.519 

.159 

-.020 

0 

0 

0 

0 

0 

0 

.039 

-.312 

.727 

.312 

-.039 

0 

0 

0 

0 

0 

0 

. 0 64 

-.516 

.933 

.516 

-.064 

0 

0 

0 

0 

0 

0 

.096 

-770 

1.142 

.770 

-.096 

0 

0 

0 

0 

0 

0 

.134 

-1.076 

1.350 

1.076 

-.134 

0 

0 

0 

0 

0 

0 

.179 

-1.432 

• 1.558 

1.432 

-.179 

0 

0 

0 

0 

0 

-.230 

1.379 

-4.138 

4.065 

.690 

0 

0 

0 

0 

0 

.862 

- 4.595 

10.338 

-13.785 

9.153 


(B 7 ) 


(B8) 


The product |~b] [^c] for use in equation (17) is found to be 


[B] M = 


.o 4 o 

.071 

.109 

.095 

.091 

.005 

.392 

-.736 

-952 

-.863 

.037 

.070 

.110 

.094 

.091 

.005 

.392 

-.736 

-952 

-.865 

-033 

.064 

.110 

.094 

.091 

.005 

.392 

-.736 

-952 

-.863 

.029 

-055 

.100 

•093 

.091 

.005 

.392 

-.736 

.952 

-.865 

.024 

.045 

.085 

.084 

.089 

.005 

.392 

-.736 

-952 

-.863 

.019 

.035 

.067 

.067 

.076 

.002 

.393 

-.736 

•952 

-.863 

.014 

.026 

.050 

.050 

.056 

-.013 

.402 

-.760 

.983 

-.891 

.014 

.026 

.050 

.050 

.056 

-.033 

-509 

-.975 

1.274 

- 1.142 

.014 

.026 

.050 

.050 

.056 

-.044 

.566 

- 1.099 

1.446 

-1.272 

.014 

.026 

.050 

.050 

.056 

-.047 

.585 

-l.l 4 l 

1-503 

- 1.310 


The dominant root of this matrix is found by iteration and is equal to 
-0.327- The normalized slope mode is given by 
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'l.OOO' 

.989 

.962 

.888 


.761, 

.575 

. 3^7 

.275 

.240 

.233j 


The dynamic pressure is given as 


1 


1 

2*( Ax) (-0.327) 


= 0.57 lb/sq. in. 


= 82 lb/sq. ft 


(BIO) 
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TABLE I 


MASS AND THICKNESS OF MODEIS 


Model 

Mass, slugs 

Thickness, in. 

31 

0. 000732 

0.016 

33 

.000928 

.020 

35 

.001446 

.032 

37 

.001875 

.o4o 

39 

.002371 

.051 

4i 

.002890 

.064 

43 

.003095 

.072 

45 

.003650 

.080 

47 

.003895 

.091 


CONFIDENTIAL 



26 


CONFIDENTIAL 


MCA RM L58E07 


TABLE II 

STIFFNESS PROPERTIES OF MODELS AND SPRINGS WITH CALCULATED 
DIVERGENCE DYNAMIC PRESSURE 


Models 

d, ft/lb 

h, ft/lb 

K e , lb/ft 

Calculated 
divergence q, 
( low-aspect- 
ratio theory) 

(a) 

Basic control surface -spring combinations 

35-55 

19.50 x 10" 5 

31.73 x 10~ 5 

19.51 

82.6 

37-58 

9-25 

14.17 

42.70 

168 

4l-6l 

2.42 

4.34 

148.1 

646 

41-67 

2.42 

1.75 

240.0 

826 

43-63 

1.71 

3.08 

208.7 

906 

43-71 

1.71 

1.17 

347.8 

1,190 - 

39-59 

4.84 

6.84 

85.7 

349 

45-65 

1.21 

3.63 

275-8 

1,222 

47-67 

.83 

2.58 

388.0 

1,762 

Modified control surface-spring combinations 

35-58 

19.50 x 10' 5 

14.17 x 10~ 5 

29.70 

103 

35-73 

19.50 

.67 

. 49.58 

l4l 

47-55 

.83 

31.73 

30.68 

652 

35-49 

19.50 

255.00 

3.64 

41.7 

31-55 

155.70 

31.73 

5.34 

16.0 

33-55 

79.80 

31.73 

8.98 

28.6 

37-55 

92.50 

31.73 

24.4 

129 

39-55 

4.93 

31.73 

27.4 

208 


divergence q calculated using piston theory is approximately 
equal to divergence q calculated by low-aspect-ratio theory multi- 
plied hy O.906M. 
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Figure 2.- Model mount used in the Langley 2-foot transonic tunnel. L-57-1783*! 
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(a) View looking toward root. L- 57 -I 782 .I 

Figure 3 . - Model mount used in the Langley 9 - "by l 8 -inch supersonic flutter tunnel. 
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(b ) View looking toward tip. L-57-1781.1 
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L- 57-1^37.1 

Figure 4.- Model mounted in 9~ by l8-inch supersonic tunnel. 



L- 57-1430.1 

Figure 5*- Model mounted in 2-foot transonic flutter tunnel. 
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L- 58-1626 

Figure 6.- Enlargements from high-speed motion picture of model 35-55 

during divergence. 
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Figure 7 .- Comparison of calculated and measured deflection modes during 

divergence . 
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Figure 8 .- Typical damaged models after supersonic testing. L-57-1438 
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Figure 9.- Variation of stiffness-altitude divergence parameter with 
Mach number for models having springs simulating symmetric mode. 
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Figure 10.- Variation, of stiffness -altitude divergence parameter with 
Mach number for models having springs simulating antisymmetric mode. 
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Effective stiffness, K e , Ib/ft 

Figure 11.- Variation of dynamic pressure at divergence with stiffness 
for model number 35 with various springs. 
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